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We introduce an adiabatic perturbation tlieory for quantum systems witli degenerate energy 
spectra. This perturbative series enables one to rigorously establish conditions for the validity of 
the adiabatic theorem of quantum mechanics for degenerate systems. The same formalism can 
be used to find non-adiabatic corrections to the non-Abelian Wilczek-Zee geometric phase. These 
corrections are relevant to assess the validity of the practical implementation of the concept of 
fractional exchange statistics. We illustrate the formalism by exactly solving a time-dependent 
problem and comparing its solution to the perturbative one. 
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Controlling the dynamics of micro and nanoscopic 
systems is at the heart of experimental and theoreti- 
cal schemes in quantum information and nanosciences. 
When, in a given physical process, there is a clear sep- 
aration of time scales into fast and slow that fact helps 
in understanding its dynamics. This is the core motiva- 
tion behind the well-known adiabatic theorem of quan- 
tum mechanics for non- degenerate systems But, how 
fast is fast and how slow is slow? It is crucial to quan- 
tify the validity of any time-dependent approximation in 
a general context since what is at stake is the practical 
dynamical implementation of concepts such as holonomic 
quantum computation or the detection of fractional (adi- 
abatic) statistics. 

A main goal of the present Letter is to asses the validity 
of the adiabatic approximation/theorem for degenerate 
spectra by introducing an Adiabatic Perturbation Theory 
[2] for Degenerate subspaces (DAPT). DAPT is a gen- 
uine perturbative expansion in the "velocity" v to drive 
the system from one particular configuration to another, 
with the perturbative corrections to the time-dependent 
Schrodinger equation (SE) given about the adiabatic ap- 
proximation (AA). Also, DAPT is not just a contribution 
to the general formalism of quantum mechanics but also a 
practical and operational formulation to test the validity 
of a broad range of concepts such as non-adiabatic frac- 
tional statistics and the extent to which AA is valid for 
a degenerate physical system, furnishing, as a bonus, the 
corrections to AA. Finally, in the development of DAPT 
we also determine corrections to the Wilczek-Zee (WZ) 
non-Abelian geometric phase Q. 

The dynamics of a quantum system with Hamiltonian 
H(t) is dictated by the SE 

i;iHvi/(,)) =H(s)|vI/(s)), (1) 

written in terms of the rescaled time variable s = vt 
(s G [0, 1]), where H = h/2TT with h denoting Planck's 
constant, T = 1/v is the relevant time scale of H(s), 
and the dot means d/ds. Consider an explicitly time- 
dependent Hamiltonian with orthonormal eigenvectors 



(s)), where g„ = 0, 1, • • • , (i„ — 1 labels states of the 
degenerate eigenspace Hn of dimension d„ and eigenen- 
ergy i?„(s), H(s)|n^'" (s)) = i5„(s)|nf" (s)), and assume 
that dn does not change during the time evolution. In 
general, Eq. ([T]) cannot be solved in closed analytic form, 
so one must resort to approximation schemes. Here, we 
are interested in developing a perturbative expansion, 
i.e., a DAPT, in the small parameter v. 

Inspired by the non-degenerate perturbative scheme 
developed in we introduce the most general ansatz 
state l^'(s)) describing a degenerate system 

ri=0 g„=0 

with n = representing the ground eigenspace, the dy- 
namical phase being given by 

W„(.S) = ^ E„is')ds' = VUJn{t), (3) 

and U"(s) a unitary matrix whose physical mean- 
ing will be elucidated below. (The un-summed in- 
dex hn is associated to the initial conditions.) By 
replacing this ansatz into ([T]) and left multiplying it 

by (m'=™(s)| one gets bn{s)Uf^^gJs) + bn{s)U'^^g^{s) + 

with w,„„(s) = w„i(s) - LUn{s), and MJ^l"g^{s) = 
(ri''" (s)|to^'" (s)), which for n^m results in 

^hZM = (n""(s)|H(s)|m3-(s))/A,„„(s), (4) 

with Am„(s) = Em{s) - En{s). 

Let us start analyzing the limiting situation v — > 0, 
defining the degenerate adiabatic approximation (DAA). 
The DAA consists in neglecting the coupling between 
different eigenspaccs "Hn but not those within a given 
eigenspace, i.e., we set 

=^n,nM^X{s) and 6„(s) = (5) 

as the degenerate adiabatic condition, implying 

f>L.„ («) + E ^z^^. (^) = and (6) 
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= E E c->"(^)[/;:'„,„ (5)&„(0)|n«"(.)). (7) 

n=0 g„=0 

In particular, if the system starts at the ground state 
|0"(0)) we have &„(0) = 6on and UJi^g^ (0) = Sh„g^ leading 
to 

\¥"Hs)) = e-^-°(-^)f/°^„(s)|0^«(s)), (8) 

where we have chosen /iq = as the initial condition 
|\I/(0)) = |0°(0)). Similar expressions hold for h„ = 
1, 2, . . . , (i„ — 1. For n = they represent the other possi- 
ble initial conditions |0i(0)), |02(0)), . . . , |0'^''-i(0)). Had 
we started in a superposition of those states the condi- 
tion UJ^^g^{0) = Sh^g„, or U"(0) — 1 using matrix nota- 
tion, must be relaxed without, nevertheless, destroying 
the unitarity of U"(s). Note that for definiteness, and 
without loss of generality, we choose /i„ = hm = 0, Vn, m, 
in the rest of this Letter. Other choices would result in 
different initial conditions for a fixed U"(0) being, never- 
theless, as good as long as we stick with it. As we will see 
below, and hinted by Eq. (g]), the adiabaticity condition 
for degenerate systems is connected to the gap between 
eigenenergies associated to each degenerate eigenspace. 

The adiabatic phase for degenerate systems has a geo- 
metric (WZ phase) and a dynamical component Q . Here 
we prove how the WZ geometric phase corresponds to 
the lowest-order approximation in our DAPT. Defining 
A^™^(s) = Eq. dSl) can be written as 

U"(s) — U"(s)A"'"(s) = (the minus sign comes from 
MJ^^g^Js) = ~{M'J^^\Js))*). By time integration of this 
expression we arrive at the WZ phase 

U"(s) = U"(0)rcxp (^J^ A""(s')ds'^ , (9) 

with T denoting a time-ordered exponential. Under a 
change of basis, i.e. a gauge transformation, U"(s) trans- 
forms unitarily after a cyclic path j^-Q ■ Note that in Ref. 
the authors assumed the initial state to correspond to 
a single eigenvector; here we have relaxed this assump- 
tion (cf. Eq. ([2])), allowing us to get a stronger condition 
for the validity of the DAA in its full generality, as we 
will see below. 

We are now in a position to state the adiabatic theorem 
for degenerate systems: "If the Hamiltonian of a system 
H(t) changes slowly during the course of time, say from 
t ~ to t ~ T, and the system is initially prepared in an 
eigenstate of H, say |n3"(0)), then it will remain in the 
instantaneous (snapshot) eigenspace Hn of H(t) during 
the interval t E [0, T]." The conditions for the validity of 
this theorem and the determination of corrections to the 
DAA when it fails are what concern us in the rest of the 
Letter. To this end, we will need to develop a DAPT. 

We want next to develop a series expansion in terms 
of the small adiabatic parameter v = 1/T, i.e., a DAPT. 



As in the non-degenerate case (2|, the choice of ansatz 
for the state |^'(s)) is crucial for its success since one 
needs to isolate the singular terms. This can be seen 
by noting that the perfect ansatz would factor out the 
dependence of |vl'(s)) on all terms of order 0{v^) and 
below. In particular, terms 0{v^^) and below are the 
problematic ones when w — >■ and should be handled 
with care. Since those terms come from e~"'^"^*), while 
the zeroth order term comes from the WZ phase, we write 
down the following ansatz 

CJO 

|*(s)) =^wP|*(P)(s)), where (10) 

p=0 

|*(p)(s)) = ^c-^""(^)B(f)(s)|n(s)) and (11) 

ri=0 

B^r^s) = ^ e^--(^)B(^),(s). (12) 

Tn=0 

Here |^(s)), j^''-^-' (s)), and |n(s)) are (i„-component col- 
umn vectors while Bi^^(s) and Bm]i(s) arc square ma- 
trices of dimensions dn x c?„, where for definiteness we 
assume dn to represent the degree of degeneracy of the 
most degenerate eigenspace. For example, |n(s))-^ = 
( |7iO(s)), |ni(s)), ••• , |n'*"-i(s)) ) . Whenever the di- 
mension of a given eigenspace m is smaller than dn, i.e., 
dm < rfrij we pad zeros to those entries and to the corre- 
sponding ones of Bmri(s)- Wc will also need the bra row 
vector, (n(5)| = ( (nO(s)|, {n\s% {n''^^-\s)\) . 

Notice that (n(s)|"^ is a column vector which combined 
to |n(s))"^ by the usual matrix multiplication rule gives 
(n(s)|-^|n(s))-^, a d„ x c?„ matrix. Each component of the 
vector |'4'(s)) has to be interpreted as corresponding to 
the evolution of an initial condition. Thus, if one starts at 
s = in the state (7„ = 1 of the ground eigenspace n = 0, 
i.e., |0^(0)), the relevant component describing its time 
evolved state is the second one, i.e., [|'i'(s))]i. The rea- 
son for adopting this convention is one of notational and 
conceptual convenience and of mathematical simplicity 
in the following calculations. We will sec that this is the 
most general ansatz accommodating an arbitrary initial 
state. 

Another important ingredient [2] for the success of 
DAPT is the setting of the right initial conditions, which 
amounts to guaranteeing the following two constraints: 

(i) The zeroth order must be the AA. In vector nota- 
tion AA reads (cf. Eq. (O), 

l*^°^(s)> = E^""""^'^^"WU"(^)|n(s)): (13) 

n=0 

implying Bi°^(s) = fo„(0)U"(s)(5„„ and b1''^(s) 
fe„(0)U"(s). Moreover, if we start at |0°(0)), implying 
fe„(0) = (5„o and U"(0) = 1, we reproduce Eq. ©. 
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(ii) For p > 1 we must have |^'^p^(0)) = 0, implying 

B(f^(0) = -E™=o,™#„B^^Uo), and thus, B^^\o) = 0. 

Having estabhshed the initial conditions one needs to 
insert Eqs. and namely, |*(s))^ = 

E.,„=oEr=o^'^e-t"-(^)|k(.rB|^i(.)^, into the SE 
(H]) in vector notation, ihv\'i{s))^ = H(s)|*(s))^, to 
obtain the following recursive relation after left multi- 
plying it with (n(s)l"^, 



fe=0 



where [M''"{s)]h 



{n<i-{s)\k^-{s)) ^ MlXis) (cf. 



]\/rnk 
^"0,1 



A,fnk 



Eq. (HI) is the following square matrix 
/ 

M'="(s) = 

(14) 

By construction, the zeroth order term is the DAA. 
The case p = and n = m, after using the initial con- 
straint (i), leads to 6„(0)(U"(s) + U"(s)M""(s)) = 0, 
thus recovering the expression for the WZ phase. 

This recursive relation is the starting point for a 
DAPT. The first non-trivial non-adiabatic contribution 
is first order in v 



•(^)6„(0)J"'""(s)U"(s)|n(s)) 



n,m— 



... ^ e-i-^.-.>„(0) -^7'°i""M |..W) 



n,m— 
m^7i 



n,m— 



■<-'^,.(.) ""''''":;'-' k.)). (15) 



with iW'""(s) = U"(s)M""(s) (U"(s))^ 

2W"™"(s') 

^nm(^ ) 



^(.)^^^d.{: 



(16) 



and 2W"™"(s) = U"(s)M"™(s)M™(s) (U"(s))^ No- 
tice that |'S''^^)(0)) = and if we have no degeneracy we 
get back to the results of Q ■ In the particular case that 
the initial state is in the ground state subspace |O(0)) 



|*(i)(s)) = ift^e-^"»("'j"""(.s)U"(s)|0(.s)) 

n=l 

.(.)iZ!!m::(f)|n(,)) 



A„o(0) 



^-..„(,) U°(.)M°»(.) 
A„o(s) 



n(.)). 



(17) 



More specifically, for an initial state |0°(0)), the first- 
order time-evolved quantum state corresponds to the first 
component l^'^^Hs)) = [l*^^Hs))]o- 

Finally, we can address the validity of the general adi- 
abatic theorem for degenerate systems. As in the non- 
degenerate case, the AA is reliable as long as the contri- 
butions coming from the first order correction are small. 
This implies the following conditions if we start at the 
ground state |0°(0)), 



hv 



n=l 



^ 1, and for all n > 1, 



hv 



[UO(s)MO"(s)] 



„nfe) [lW°"(0)U"(5)]o3„ 



A„o(s) 



A„o(0) 



«1, 



reducing to the ones in 0] when there is no degeneracy. 

To illustrate our DAPT formalism, consider a four- 
level quantum system subjected to a rotating classical 
magnetic field B(t) = Br{t) of constant magnitude B. In 
spherical coordinates r{t) = {s'md cos (f{t), s'm9 smif{t), 
cos 9), with < 9 < TT and < < 27r the polar and az- 
imuthal angles^respectively. The Hamiltonian describing 
the system is Q H(t) = f 6r(t) • T, with 6 > propor- 
tional to B and F = {Tx, Ty, P^) Dirac matrices satisfy- 
ing the Clifford algebra {P^, Tj} = 2(S.yl4, j = x, y, z. In 
terms of Pauli matrices aj , we choose the particular rep- 
resentation Pj = cFx®cFj, so that [Pi,Pj] = 2ieijfcnfc, with 
Ilj = I2 ® (Jj , and eijk the Levi-Civita symbol. Then, the 
snapshot eigenvectors in the basis where II^ is diagonal 
are (n = 0, 1 and a = e"^*^*-' sin0, (3 = cos 6*) 



1 



\n\t))^^{a*\n)-P\n) 



1)"IU)), (18) 
(19) 



\n\t))^^{m)+c^\n)-{-mm. 



with two-fold degenerate, time-independent, snapshot 
eigenvalues Eq = —{h/2)b and Ei = (h/2)h. 

We can exactly solve the time-dependent SE ([1]) when 
Lp{t) = wt, w > representing the frequency of the ro- 
tating magnetic field, by employing techniques similar to 
those developed for the single spin- 1/2 problem iS-jl^l. 
Assuming that at t = the initial state is |0°(0)), the 
resulting time-dependent solution, expressed in terms of 
the snapshot eigenvectors of H(t), is 



e' 2 



1 + cos [ 



■A-it) 



1 — cos ( 



■A+it) 



|0°(0) 



„ A+(t) - A_(t), 

e^'— sin^ +^ ' ^10^(0) 



e 2 sm 



.^,^B,it) + B^it) 



sm ( 



2 

1 + cos ( 



{t)) 



1 — cos 9 , . 
■B-it) B+it) 



4 



biw cos 9 



O4 



sin(rj±t/2), B±(t) 



2 + 62 ± 2 w 6 cose*. 



with A±{t) = cos (r2±t/2) + i 

i sin(rj±t/2), and = w 

In order to develop a DAPT expansion one needs the 
expression for the WZ phase. For our problem it can be 
determined exactly (n = 0,1) 



U"(t) 



Z2 ^1 



, (20) 



where Zi = e™*/^ [cos cos6i) - icos6'sin cos^)] 
and Z2 = ie™*/2 gin^sin cos6'). In this way, the ze- 
roth order term, i.e., the DAA reads 

|vl/(0)(t)) = e-^-«W(f/°o(t)|0"(t))+t/„",(t)|Oi(t))) 

= e'^(zi|0°(i))-^2|0^W)) (21) 

after replacing ujo{s) = —bs/2 ~ — &wt/2, which is the 
same obtained by expanding the exact solution |4'(t)) up 
to zeroth order. The first non-trivial correction is first 
order in DAPT, where the small parameter is v = w. It 
is given by 

|^r(i)(t)) = i^sin2e'|*^"Vt)) +i^sin(6t/2)sin6' 
ibv bv 

X (zi sin 6* + Z2 cos 6*) 1 1" (t )) + ^ [cos(5t/2)z2 

bv 



isin(6i/2) cos 6'(zi* 



Z9 cos t 



ll'(^)>, 



(22) 



and it exactly coincides with the first order expansion of 
the exact solution |4'(t)). Note that in working out the 
expansion of the exact solution one must be aware that 
terms w"^^t are actually order w" since t (x 1/v = 1/w. 

We should emphasize that for most time dependent 
Hamiltonians one is not able the get the exact solution. 
However, our DAPT is general enough to provide a sys- 
tematic expansion in powers of v about the DAA for 
any time dependent degenerate Hamiltonian. Further- 
more, DAPT can be used to get corrections to the non- 
abelian WZ phase. In our example, if the first order 
correction to the DAA becomes relevant, the system's 
state can be written as |*(s))Ar = ^P(s)|$(°)(s)) + 
^g-ifat/2]v(s)(co(s)|lO(s)) + ci{s)\l\s))) + 0{v^), where 
N{s) normalizes the state, P (s) = |(*(°)(s)|*(s))Arp, 
|$(0)(s)) = N{s)[l + t;/(s)]//P(i)|vI/(o)(s)), and /(s) = 
is sin^ 0/(46). The quantities co(s) and ci(s) are not 
needed and can be obtained via Eq. ((22|). P{s) gives the 
probability of measuring the first order corrected state 
|^(s))7V in tfic adiabatic state |\I'''''(s)), i.e., the chances 
of being at ground eigcnspace. The remaining terms in 
|$(°ns)) determines how |0"(s)) and \0^{s)) are now su- 
perposed within the ground eigcnspace. |$^'^)(s)) is what 
is actually seen by any experiment trying to get the WZ 
phase if the first order is relevant. Keeping terms up 
to first order we get |$(°)(s)) = [1 -I- w/(s)]|1'(°)(s)> = 
e'T {xi |00(s)) + X2 |Oi(s))), where xi = {I + vf{s))zi 



and X2 = —(1 + vf{s))z2 (cf. Eq. (|2T|) V We can repeat 
the previous argument but with the system starting at 
|0^(0)), which gives the two remaining terms for the cor- 
rected WZ phase (cf. the second row of Eq. (f20| with 
n = 0). Putting everything together and noting that v = 
w we get instead of the WZ phase U^°^ [t) the following 
non-abelian phase V(°)(t) = U(")(t) + i«-"t^m^ (t). 
This is what is experimentally seen if the first order cor- 
rection is relevant. Note that V^'^^ (i) is unitary up to first 
order and that for w — )• we recover continuously the WZ 
phase. Also, the previous procedure applied to the non- 
degenerate system of [H gives the same correction to the 
Berry phase therein computed by another method. 

In conclusion, we have shown a genuine adiabatic per- 
turbation theory for Hamiltonians with degenerate spec- 
trum (DAPT) as a series expansion about the degenerate 
adiabatic approximation (DAA), in terms of the "veloc- 
ity" V by which the system is driven away from its initial 
configuration. DAPT allowed us to get a general formu- 
lation to the adiabatic theorem for degenerate systems as 
well as corrections to the non-abelian Wilczek-Zee phase 
when DAA no longer holds. Finally, the key ingredients 
in the construction of DAPT, differentiating it from all 
standard perturbation theories and other approaches @ , 
are three fold: (a) the rescaling of the real time t to 
s = vt, allowing a consistent perturbative expansion to 
all orders in terms of v\ (b) the correct setting of the ini- 
tial condition at t = 0, i.e., one should recover the zeroth 
order not only when w — > but also at < = 0; and (c) 
the powerful vectorial ansatz, Eqs. (fT0 |) - (fT2|) . that fac- 
tored out singular terms and permitted the construction 
of simple matrix recursive relations. 
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